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ABSTRACT 



We use 3D cosmological simulations with radiative transfer to study the formation 
and evolution of the first galaxies in a AC DM cosmology. The simulations include 
continuum radiative transfer using the "Optically Thin Variable Eddington Tensor" 
(OTVET) approximation and line-radiative transfer in the H2 Lyman- Werner bands of 
the UV background radiation. Chemical and thermal processes are treated in detail, 
particularly the ones relevant for H2 formation and destruction. 

We find that the first luminous objects ( "small-halos" ) are characterized by bursting 
star formation (SF) that is self-regulated by a feedback process acting on cosmological 
instead of galactic scales. The global star formation history is regulated by the mean 
number of ionizing photons that escape from each source, ejjvifesc)- It is almost in- 
dependent of the assumed star formation efficiency parameter, e*, and the intensity 
of the dissociating background. The main feedback process that regulates the SF is 
the re- formation of H2 in front of Hii regions and inside relic Hii regions. The Hii 
regions remain confined inside filaments, maximizing the production of II2 in over dense 
regions through cyclic destruction/reformation of II2. If ejjvifesc) > 10^ ''/e* the SF is 
self-regulated, photo-evaporation of "small-halo" objects dominate the metal pollution 
of the low density IGM, and the mass of produced metals depends only on {fesc)- If 
euvifesc) ~ 10~^/e*, positive feedback dominates, and "small-halo" objects constitute 
the bulk of the mass in stars and metals at least until redshift z ~ 10. "Small-halo" 
objects cannot reionize the universe because the feedback mechanism confines the H 11 
regions inside the large scale structure filaments. In contrast to massive objects ("large 
halos"), which can reionize voids, "small-halo" objects partially ionize only the dense 
filaments while leaving the voids mostly neutral. 

Subject headings: Cosmology: early universe — cosmology: theory — galaxies: dwarf — 
galaxies: evolution — galaxies: formation — galaxies: high-redshift — intergalactic medium — 
methods: numerical 



DRAFT: February 1, 2008 



2 



1. Introduction 

In cold dark matter (CDM) cosmologies, "small-halo" galaxies are believed to be the first 
luminous objects formed in the universe. Defined as small mass protogalaxies, "small-halo" objects 
have virial temperature Tyir < 10^ K and rely on H2 line cooling to form stars. The first generation 
of stars is necessarily metal-free, as all elements heavier than Li are produced in the cores of stars 
or by supernova (SN) explosions. In the literature Population III is often used to refer to both 
metal-free stars and to protogalaxies with T^r < 10^ K. In order to avoid this confusion, we use 
the term "small-halo" objects instead of the widely used "Population III objects" throughout our 
paper. For instance, the stars in a "small-halo" object are not necessarily Population III because 
the interstellar medium (ISM) could be polluted by metals rather quickly. In this paper, we use 
synthetic stellar energy distributions (SEDs) calculated for metal-free stars, which for brevity we 
call Population III SEDs, and SEDs typical of low-metallicity stars, which we call Population II 
SEDs. Although the SED does not change the star formation (hereafter SF) history, we find that 
the H2 abundance in the intergalactic medium (IGM) depends on it. 

In this paper, the second in a series on the formation and evolution of the first galaxies, we 
focus on the mechanisms responsible for the the self-regulation of star formation in "small-halo" 
objects. In Ricotti, Gnedin, &; Shull (2002) (Paper I) we discussed the numerical methods and 
the physics included in the simulations used in this paper. From resolution studies we found that 
our largest simulations seem to be close to numerical convergence. These models are the first 3D 
cosmological simulations to include the physics necessary to study the formation of "small-halo" 
objects and their radiative feedback self-consistently. Mechanical feedback from SN explosions and 
stellar winds is not included, as it is the subject of a separate paper. The main new ingredients are 
the inclusion of continuum radiative transfer and line radiative transfer in the Lyman- Werner bands. 
A number of potentially relevant physical processes are included as well: secondary ionizations of H 
and He, detailed H2 chemistry and cooling processes, heating by Lya resonant scattering, H and He 
recombination lines, metal production, and radiative cooling. The SED of the sources is consistent 
with the choice of the escape fraction of ionizing radiation, {fesc)- The results of our simulations 
are considerably different from previous models of the formation of "small-halo" objects. These 
previous works were semi-analytic treatments (Haiman, Rees, & Loeb 1997; Haiman, Abel, & Rees 
2000; Ciardi et al. 2000) or 3D simulations (Machacek, Bryan, & Abel 2001) without radiative 
transfer. In contrast to previous studies, we find that the formation of "small-halo" objects is not 
suppressed by the H2 dissociating background. The reason for this difference is that these earlier 
studies did not include a primary mechanism of positive feedback for the formation of "small-halo" 
objects. 

Prior to our work, it was widely thought that the main feedback mechanism that regulates 
"small-halo" formation is the build-up of the H2 dissociating background. Indeed, the dissociating 
background can suppress or delay the formation of "small-halo" objects if H2 is not reformed 
efficiently. Because of these results, "small-halo" objects were believed to be unimportant for 
subsequent cosmic evolution, metal enrichment of the IGM, and reionization. However, our recent 
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work (Ricotti, Gnedin, & Shull 2001) pointed out a new radiative feedback mechanism that turns 
out to provide the dominant regulation for the formation of "small-halo" objects. A large amount of 
H2 is naturally reformed in shells (PFRs: positive feedback regions) in front of Hll regions or inside 
recombined fossil H 11 regions. The H2 formation rate is proportional to the density squared of the 
gas and is very efficient in the filaments and inside galaxies. The bursting mode of SF observed in 
simulations produces fossil H II regions and therefore allows H2 to be continuously reformed. 

Contrary to the initial ideas in Ricotti et al. (2001), the volume filling factor of positive 
feedback regions (PFRs) has to remain small to maximize the positive feedback. Therefore, we do 
not require ionizing escape fraction {fesc) ~ 1? to maximize the IGM volume occupied by PFRs with 
respect to the volume occupied by the dissociating spheres (or with respect to the intensity of the 
dissociating background). Instead, "small-halo" formation is possible if {fesc) < 1; inside the dense 
filaments, the positive feedback of H2 reformation always dominates the negative feedback of the 
dissociating background and the of local dissociating radiation. This, together with the clustering 
of dark matter (DM) halos in the over dense regions, maintains a sufficiently high H2 abundance 
inside the filaments to allow "small-halo" objects to cool and form stars. On the other hand, if 
the star formation rate (SFR) is too high (and (fesc) ^ 1), the filaments become highly ionized 
by the Stromgrcn spheres and the H2 is destroyed by direct ionization. This produce a local and 
temporary halt of SF until the Stromgren spheres recombine. Shells of H2 then reform inside the 
H II regions and a new burst of SF occurs. 

In Ricotti et al. (2001) we emphasized the importance of PFRs in reducing the dissociating 
background. In some simulations, the production of H2 is high enough to reduce the dissociating 
background intensity by a factor of 10. But this does not affect the SF history. In the same paper, 
we pointed out that enhanced galaxy formation is possible if it happens inside a PFR. Indeed, 
in this work we find that this is the most important mechanism that regulates the formation of 
"small-halo" objects. Another positive feedback mechanism proposed by Fcrrara (1998) involves 
local production of H2 in shells produced by SN explosions. In the simulations presented in this 
paper, we do not include feedback from SN explosions, although we plan to address this problem 
in a subsequent paper. 

In Ricotti & Shull (2000) we studied (fesc) for "small-halo" objects, assuming a spherical 
halo geometry. The high merger rate in the early universe should favor the formation of spheroidal 
galaxies instead of disks. We have found that {fesc) is small at high redshift, {fesc) e^p[—{l+Zyir)]. 
This supports the theory that "small-halo" objects formed copiously in the early universe, and 
their existence could be directly or indirectly observed today. Metals in the Lya forest and dwarf 
spheroidal galaxies in the Local Group are two examples of observations that can test our models. 
In Paper III, currently in preparation, we will study the properties of "small-halo" galaxies and 
try to make a connection with available observational data on the dwarf spheroidal galaxies in the 
Local Group. 

The paper is organized in the following manner. In § 2 we briefly review the physics included 
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in the code and the free parameters. This section serves as a quick reference of arguments treated 
extensively in Paper I . The results of the simulations are shown and discussed in § 3. In § 4 we 
provide a summary and final comments. 

2. The Code 

The simulations were performed with the "Softened Lagrangian Hydrodynamics" (SLH-P^M) 
code described in detail in Gnedin (1995, 1996), and Gnedin &; Bertschinger (1996). The simulation 
evolves coUisionless DM particles, gas particles, "star-particles" formed using the Schmidt law in 
resolution elements that sink below the numerical resolution of the code, and radiation, whose 
transfer is treated self-consistently with the OTVET approximation of Gnedin & Abel (2001). We 
also include line radiative transfer in the H2 Lyman- Werner bands of the background radiation, 
secondary ionization of H and He, heating by Lya scattering, detailed H2 chemistry and cooling, 
and self-consistent SED of the sources (Paper I). 

We adopt a ACDM cosmological model with parameters: r^o = 0.3, = 0.7, h = 0.7 and 
Qj, = 0.04. The initial spectrum of perturbations has as = 0.91 and n = 1. All simulations start at 
z = ICQ and finish at z ^9. We use box sizes L(,oa; = 0.5, 1 and 2 comoving h^^ Mpc and grids with 
^box = 256^, 128^ and 64^ cells. We achieve the maximum mass resolution of Mdm = 4.93 x 10^ 

Mq and spatial resolution of 156 comoving pc in our biggest run. We fully resolve the SF 
in objects within the mass range 5 x 10^ Mq ^ Mdm ~ 10^ Mq. 

In Paper I we discussed extensively the details of the code and the physics included in the 
simulation. We also studied the numerical convergence of the simulations that is especially crucial 
in the study of the first objects. High mass resolution is needed because the objects that we 
want to resolve have small masses (10^ Mq<Mdm^ 10^ Mq). Moreover, the box size has to 
be large enough in order to include at least a few of the rare first objects. The first "small-halo" 
objects should form at z ^ 30 from 3a density perturbations; the first "large-halo" objects (with 
Mdm ^ 5 X 10^ h"^ M0) should form at z ^ 20, also from 3a perturbations. 

The reader interested in numerical issues or in the details of the physics included in the 
simulation should refer to Paper I. Below, we summarize the meaning of the four free parameters 
in the simulations. 

• e*: Star formation efficiency in the Schmidt law (dp^/dt = e^^pg/t^,, where /3* and pg are 
the stellar and gas density, respectively, and is the maximum between the dynamical and 
cooling time). 

• euv- Energy in ionizing photons per rest mass energy of H atoms [mnc^) transformed into 
stars. This parameter depends on the IMF and stellar metallicity. 

• if esc)'- escape fraction of ionizing photons from the resolution element. 
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• g,y: Normalized stellar energy distribution (SED). We use Population III and Population II 
SEDs with Salpeter (1 < < 100 Mq) initial mass function (IMF), and we modify gi, 
according to the value of (fesc)-^ 

3. Results 

In this section, we discuss the main physical processes that regulate the formation and evolution 
of "small-halo" objects. We ran a large set (~ 20) of simulations in order to explore the dependence 
of the results on free parameters in the simulation. We ran the simulations on the Origins2000 
supercomputer at the NCSA in Urbana-Champaign, IL. The typical clock time to run a 64^, 128^, 
and 256^ simulation from z = 100 to z = 10 is about 50, 1000, and > 23, 000 hours. The total 
computational time used to run the simulations presented in this work is about 44, 000 hours. In 
Table 1, as a quick reference, we list the simulations with radiative transfer discussed in this section. 
Each simulation is named, 64L05p2noLW-|-l for example, using the following convention: 

^L05[p2j [noIM] [+lj 

Nbox / i \ ^ Oclog(e[/y(/esc)) 

Lbox Pop comment 

where Nj^^^ is the number of cells in the box, Li,ox is the comoving size of the box in h^^ Mpc, "p2" 
refers to Population II (metallicity Z = 0.04 Zq) and "p3" to Population III (metal-free) stars. The 
brackets indicate the optional part of the name. The comment "noLW" indicates that the sources 
emit no Lyman- Werner H2 photons, "noRAD" is used for simulations without radiative transfer, 
and "noC" indicates that radiative transfer and H2 cooling are not included in the simulation. The 
second optional part of the name is log[(e[/y/4-7r)(/esc)/l-6 x 10^^]. Values of zero are omitted from 
the name. We will attempt to make clear the parameters of each simulation in the text. Therefore, 
a reader does not need to memorize the notation or repeatedly refer to Table 1. 



small value of (/esc) increases the flux jump at the Lyman limit and makes the SED harder. 
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Tabic 1. List of siniTilations with radiative transfer. 



RUN Nbo^ Lbox Mass Res. {U^c) {^){fesc) e* (^)(/e.c)e* 

{h-^ Mpc) [h-^ M0 ) 



64L05noRAD 


64 


0.5 


3.94 


X 


10*^ 


10 


- 


1 









0.2 







64L05noC 


64 


0.5 


3.94 


X 


10*^ 


10 




1 









0.2 







o4LiUopznoLW -2 


64 


U.O 


3.94 


X 


1 n4 
lU 


10 


TT 

11 


1 


l.D 


X lU 


-7 


U.Z 


3.2 


X 10-^ 




64 


n 5 


3 Q4 






10 


II 




J- .LI 


X 10" 


-5 


n 2 


3.2 


X 10-6 




64 


n 5 


3 Q4 




J-\J 


10 


n 




1 fi 


X 10" 


-4 


02 


3.2 


X 10-6 




64 


5 


3 94 




10^ 


1 


n 




1 fi 

-L . \J 


X 10" 


-7 


n 2 


3.2 


X 10-s 


64L05p2-l 


64 


0.5 


3.94 


X 


10^ 


10 


II 


1 


1.6 


X 10" 


-6 


0.2 


3.2 


X lO-'^ 


64L05p2 


64 


0.5 


3.94 


X 


10^ 


10 


II 


1 


1.6 


X 10" 


-5 


0.2 


3.2 


X 10-6 


64L05p3 


64 


0.5 


3.94 


X 


10^ 


10 


III 


1 


1.6 


X 10" 


-5 


0.2 


3.2 


X 10-6 


64L05p3b 


64 


0.5 


3.94 


X 


10^ 


10 


III 


1 


1.6 


X 10" 


-5 


0.02 


3.2 


X 10-'^ 


64L05p3c 


64 


0.5 


3.94 


X 


10^ 


10 


III 


1 


1.6 


X 10" 


-5 


0.002 


3.2 


X 10-s 


64L05p2-lf 


64 


0.5 


3.94 


X 


10^ 


10 


IP 


0.1 


1.6 


X 10" 


-5 


0.2 


3.2 


X lO-'^ 


64L05p2-2f 


64 


0.5 


3.94 


X 


10^ 


10 


IP 


10-2 


1.6 


X 10" 


-5 


0.2 


3.2 


X 10-s 


64L05p3-2f 


64 


0.5 


3.94 


X 


10^ 


10 


iii*^ 


10-2 


1.6 


X 10" 


-5 


0.2 


3.2 


X 10-^ 


64L05p2-2fa 


64 


0.5 


3.94 


X 


10^ 


10 


IP 


10-2 


1.6 


X 10" 


-5 


0.05 


8 X 10-9 


64L05p3-2fa 


64 


0.5 


3.94 


X 


10^ 


10 


iii^ 


10-2 


1.6 


X 10" 


-5 


0.05 


8 X 10-^ 


64L05p2-3 


64 


0.5 


3.94 


X 


10^ 


10 


IP 


10-3 


1.6 


X 10" 


-5 


0.2 


3.2 


X 10-9 


64L05p2-5 


64 


0.5 


3.94 


X 


10^ 


10 


IP 


10-5 


1.6 


X 10" 


-5 


0.2 


3.2 


X 10-11 


64L05p3b-3n 


64 


0.5 


3.94 


X 


10^ 


10 


III 


1 


1.6 


X 10" 


-8 


0.02 


3.2 


X 10-10 


64L05p3b-3^ 


64 


0.5 


3.94 


X 


10^ 


10 


iii^ 


10-3 


1.6 


X 10" 


-8 


0.02 


3.2 X 10-1° 


64LlnoRAD 


64 


1.0 


3.15 


X 


10^ 


10 




1 









0.2 







64LlnoC 


64 


1.0 


3.15 


X 


10^ 


10 




1 









0.2 







64Llp2 


64 


1.0 


3.15 


X 


10^ 


10 


II 


1 


1.6 


X 10" 


-5 


0.2 


3.2 X 10-6 
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In Figure 1 (left) we show the comoving SFR (Mq yr-i Mpc-^) for the 64L05noRAD, 64L05noC, 
64L05p2, 128L05noRAD, and 128L05p2 runs. These simulations have L^ax = 0.5 Mpc and 
= 0.2. The two solid lines show the SFR in the 128^ box (thick line) and in the 64'^ box (thin 
line) without radiative transfer. In the higher mass resolution simulation, the SFR is larger since 
more small mass objects are formed. The short-dashed line shows the 64^ simulation without radia- 
tive transfer and H2 cooling. By definition, this simulation does not form any "small-halo" objects. 
Comparing the short-dashed line with the solid lines, we see that "small-halo" objects dominate the 
SFR down to 2; ~ 10 if we do not include radiative feedback. The two long-dashed lines show the 
SFR including radiative transfer (using the Population II SED and (e[/y/47r)(/esc) = 1-6 x 10~^). 
Again, the thick line is for the 128^ box and the thin line is the 64^ box. From Figure 1 it is clear 
that SF is bursting and is suppressed with respect to the simulations without radiative transfer. 




z+1 z+1 

Fig. 1.— (Left) SFR in 0.5 Mpc box simulations: 64L05noRAD, 64L05noC, 64L05p2, 
128L05noRAD, and 128L05p2. The solid lines show simulations without radiative transfer and the 
long-dashed lines including radiative transfer. The short-dashed line shows a simulation without 
either radiative transfer or H2 cooling. (Right) The same as on the left figure but for 1 Mpc box 
simulations: 64LlnoRAD, 64LlnoC, 64Llp2, 128LlnoRAD, and 128Llp2. The 64L1 simulations 
have insufficient mass resolution to form "small-halo" objects. 

Figure 1 (right) shows the analogous simulations to Figure 1 (left) but for box size L^ox = 1 
Mpc. The three thin lines show the SFR for the 64^ box: without radiative transfer (solid), 
without radiative transfer and H2 cooling (dashed), and with radiative transfer (long-dashed). 
The lines are almost indistinguishable, showing that we are forming only "large-halo" objects and 
that radiative feedback has no effect on their SFR. The DM mass resolution of these simulations, 
Mdm = 3.15 X 10^ h^^ Mq, is not sufficient to fully resolve the first "small-halo" objects, with 
typical masses Mdm ~ 10^ - 10^ h-^ Mq. In Paper I we showed that we need to resolve each 
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Table 1 — Continued 



RUN 


-^box 


■^box 

{h-^ Mpc) 


Mass Res. 
{h-^ Mo ) 




9u 


if esc) 




e* 


{ 4Tr) if esc) ^* 


128L05noRAD 


128 


0.5 


4.93 X 10^ 


10 




1 







0.2 





128L05p2 


128 


0.5 


4.93 X 10^ 


10 


II 


1 


1.6 X 10" 


-5 


0.2 


3.2 X 10-^ 


128LlnoRAD 


128 


1.0 


3.94 X 10^ 


10 




1 







0.2 





128Llp2 


128 


1.0 


3.94 X 10^ 


10 


II 


1 


1.6 X 10" 


-5 


0.2 


3.2 X 10-6 


128Llp2-2^ 


128 


1.0 


3.94 X 10^ 


16 


IP 


10-2 


1.1 X 10" 


-7 


0.05 


5.5 X 10-9 


256LlnoRAD 


256 


1.0 


4.93 X 10^ 


10 




1 







0.2 





256Llp3^ 


256 


1.0 


4.93 X 10^ 


25 




0.1 


2.5 X 10" 


-6 


0.1 


2.5 X 10"^ 



Note. — Parameter description. Numerical parameters: N^^^ is the number of grid cells, Lbox is the 
box size in comoving Mpc, is the parameter that regulates the maximum deformation of the 
Lagrangian mesh: the spatial resolution is ^ Li^y,./ {Ni),yj.B^,). Physical parameters: is the normalized 
SED (II=Population II and III=Population III), e* is the star formation efficiency, euv is the ratio of energy 
density of the ionizing radiation field to the gas rest-mass energy density converted into stars (depends on 
the IMF), and {fesc) is the escape fraction of ionizing photons from the resolution element. 

^Qj^ is modified assuming ag = iVneiZ-^Hi = 0.1 and ai = Nueu/N-ni = 10 where Ni is the column density 
of the species/ion i (see Paper I). 

^g^ is modified assuming oq = 0.01, ai = 10. 

'^Secondary ionizations included. 
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object with at least 100 DM particles. The two thick lines show the SFR for the 128^ boxes 
without radiative transfer (solid line) and including radiative transfer (long-dashed line). For the 
simulation with radiative transfer we use the Population II SED and (ec/y/47r)(/esc) = 1-6 x 10~^. 
This simulation has mass resolution Mdm = 3.94 x lO'* M©, sufficient to resolve "small-halo" 
objects, and box size large enough to include the first "large-halo" objects that form at z ~ 20. It 
is shown clearly that the star formation is bursting at z ^ 20 when it is dominated by "small-halo" 
objects. In this simulation, "large-halo" objects dominate the SFR as soon as they form, producing 
the continuous star formation mode observed at 2; ^ 20. 

In summary, simulations with mass resolution Mdm < 3.15 x 10^ Mq (such as the 64L1 
runs) do not form "small-halo" objects and therefore are not well-suited for this study. On the other 
hand, if the box is much smaller than 1 Mpc (comoving), the formation of the first "large-halo" 
objects is delayed considerably. If L^ox > 1 Mpc, the first rare massive objects form at 2; ~ 20. 
From the results of other simulations, not shown here, we have verified that the global SFR does 
not change if we increase the box size from L^ox = 1 Mpc to L^ox = 8 h"^ Mpc, keeping 
the mass resolution constant (for redshifts z > A). The bursting SF mode of "small-halo" objects 
is not synchronized throughout the whole universe. Therefore, the strong oscillations of the SFR 
observed in some of our simulations are an artifact of the finite (actually quite small) size of the 
simulation box. 

In Figure 2 we show a sequence of four 3D-views of the cube for the 128Llp2 run at z = 
21.2, 17.2, 15.7, 13.3. The rendering of the volume is obtained by assigning quadratic opacity to 
the logarithm of the gas density and linear opacity to the xhu fraction. The colors show the xhh 
fraction. The cosmological H II regions expand quickly in the IGM, but when their size becomes on 
the order of the filamentary structure the expansion stops (in comoving coordinates) . The filaments 
become partially ionized, but the voids remain neutral and reionization cannot occur. At redshift 
z ~ 20 the first "large-halo" objects form. For "large-halo" objects, there is no feedback process 
that can stop the H 11 region from expanding into the low density IGM. Indeed, these H 11 regions 
never stop expanding and eventually reionize the universe. Figure 3 is analogous to the previous 
figure, but here we show the H2 abundance for the 64L05p3 run. The H2 is quickly destroyed in 
the low density IGM by the build-up of the dissociating background, but new II2 is continuously 
reformed in the filaments. Therefore, in the dense regions where galaxy formation occurs, positive 
feedback dominates over the negative feedback of the dissociating background. In the following 
sections, we will show that the main mechanisms that regulate galaxy formation of "small-halo" 
objects are the two positive feedback processes found by Ricotti et al. (2001). In that paper, we 
discussed the importance of H2 shells (PFRs: positive feedback regions) that form just in front 
each Stromgren sphere and inside relic (recombining) Hll regions. 

Since the 128^ and 256^ cubes are computationally intensive, in the following three sections 
(§ 3.1-§ 3.3) we will try to understand the feedback mechanism that regulates the SF of "small- 
halo" objects using 64'^ cubes with L^ox = 0.5 h^^ Mpc. We are aware, though, that after z ~ 20 
we are not properly including the formation of the first "large-halo" objects. In § 3.4 we present 
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Fig. 2. — Rendering of the 128Llp2 run at z = 21.2, 17.2, 15.7, 13.3. The colors code the logarithm 
of the H II abundance. "Small-halo" objects produce H ii regions that remain confined in the denser 
filaments, creating structures that resemble a pearl necklace. At z = 13.3 the first "large-halo" 
object produces a large Hii region (big blob on the left of the image) that continues to expand and 
eventually reionizes the IGM around it. 




Fig. 3. — Rendering of the 64L05p3 run at z 
of the H2 abundance. 



= 15.7, 13.3, 12.3, 10.2. The colors code the logarithm 
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the results of our largest simulations that include detailed physics, and realistic values of the free 
parameters. 

3.1. Negative and Positive Feedback 

The background intensity in the Lyman- Werner bands determines the redshift at which H2 is 
destroyed in the low density IGM. In Figure 4 (top) we show the mean mass- and volume- weighted 
molecular abundances (x^/j)) {xh+) i^H-) ^ ^ function of redshift for the 64L05p2noLW, 
64L05p2 and 64L05p3 runs. These runs have the same parameters (e* = 0.2, (euv/^T^) = 1-6 x 10~^, 
and (fesc) = 1) except for a different flux in the Lyman- Werner bands^. The solid line shows the 
64L05p2noLW run for which the stars do not emit Lyman- Werner photons. The short-dashed line 
shows the 64L05p2 run that has a Population II SED, and the long-dashed line shows the 64L05p3 
run that has a Population III SED (zero metallicity) . The Population II SED emits about 100 times 
more photons (per Mq of SF) in the Lyman- Werner bands than the Population III SED. Indeed, 
in the 64L05p2 run, the dissociating background reduces the H2 abundance in the low-density IGM 
to a value 1/100 of the H2 abundance in the 64L05p3 run. In the 64L05p2noLW the H2 is not 
destroyed. In Figure 4 (bottom) we show the SFR as a function of redshift for the same simulations. 
Surprisingly, the SFR does not depend appreciably on the intensity of the dissociating background. 
It is evident that the destruction of H2 in the low-density IGM does not affect the global SFR. The 
only difference in SF history is the run with a Population II SED, which has the higher background 
in the Lyman- Werner bands. Here, the SFR decreases more than in the other two runs before 
a new burst occurs. However, on average, also taking into account numerical errors, the SFR is 
indistinguishable in these three runs. 

When the background in the Lyman- Werner bands (averaged over 1000-1100 A) builds up 
to Jlw ~ 10~^^ erg cm~^ s~'^ Hz"^ sr"-*^ (at z ~ 20), it starts to dissociate the H2 in the IGM. 
It is evident (Fig. 4) that H2 is continuously destroyed and reformed. The H2 abundance peaks 
are slightly delayed with respect to the SFR minima. Since SF is self-regulated, the dissociating 
background intensity, after a rapid build-up phase, reaches an almost constant value. Depending 
on the SED of the sources and (fesc): the equilibrium value of the background can be small (e.g.. 
Population HI and high (fesc)) or large (e.g.. Population II and small (/esc))- Consequently, the 
H2 abundance in the IGM can assume either a large or small quasi-constant value during the 
self-regulated SF phase. 

In the Population HI SED run (long-dashed line), the molecular abundance in the IGM is 100 
times higher than in the Population II run, and the mass-weighted abundance becomes {xh2)m ~ 
2 X 10~^ at z ~ 12. This means that the positive feedback produced essentially the same mass of 



^The SED, g^, of ionizing radiation is different for Population III and Population II, but since euv is fixed, the 
Hi ionization rate is the same. 
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Fig. 4. — (Top) Species abundances for the 64L05p2noLW (solid line), 64L05p2 (short-dashed 
line), and 64L05p3 (long-dashed line) runs. These runs have the same parameters (e* = 0.2, 
(ef/y/47r) = 1.6 x 10^^ and {fesc) = 1) except for a different flux in the Lyman- Werner bands. The 
two top panels show the volume (left) and mass (right) weighted H2 abundance in the box. The 
middle and bottom panels show the analogous quantities for and H~ respectively. (Bottom) 
The comoving SFR for the same three runs. Despite the fact that the dissociating background flux 
in these three simulations differs by several orders of magnitude, the SFRs are almost identical. 
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H2 destroyed by negative feedback. 

The opacity of the IGM in the Lyman- Werner bands is proportional to the mean H2 number 
density. In Figure 5 we show the background specific intensity in 10~^^ erg cm~^ s"^ Hz~^ sr~^ 
units, J2i('^)) SlS a function of frequency for the 64L05p3. Wc show J2i{j^) at redshift z = 12 
when {xh2)m ~ 2 x 10^^. In the upper panel, it is interesting to note the importance of the 
spectral features caused by Hi , Hei , and Hell Lya emission lines. The lower panel shows a 
zoom at the frequencies of the Lyman- Werner bands; the upper (lower) line shows the intensity of 
the background without (with) line radiative transfer. The H2 and resonant H I Lyman series line 
opacities reduce the background intensity by about one order of magnitude. 

In Figure 6 we show a time sequence of two slices through the most massive object for the 
64L05p2 simulation at z = 19.4 and 18.5. This simulation has e* = 0.2, (e(7y/47r) = 1.6 x 10~^, 
ifesc) = 1) Eind sources with a Population II SED. Each one of the two panels shows logXHi, 
log.Tj:/2 5 log(l + ^) ^-iid logT, where 6 = (p — po)/po is the baryon overdcnsity with respect to 
the mean IGM density po. At z = 19.4, the H2 has its relic abundance everywhere in the IGM 
except in the dissociation spheres around the first objects, where it is destroyed. At 2; = 18.5, the 
dissociation spheres are still visible, but the UV background starts to dissociate H2 everywhere in 
the IGM except in the filaments. Finally, at z ^ 18.5, before the dissociation spheres overlap, the 
background has destroyed all the relic H2. The H2 is still present in the filaments where the gas is 
partially ionized by stars, and positive feedback dominates. In the analogous simulation, where the 
sources have a Population III SED, the dissociation spheres never appear around the source. The 
dissociating background destroys the H2 in the IGM before the dissociation spheres grow larger than 
the PFRs. Finally, if (fesc) < 1 and we use a Population II SED for the sources, the dissociation 
spheres around the sources almost overlap before the background dominates the H2 dissociation 
rate. If the dissociating radiation emitted by each source, Slwi is large, the dissociating background 
intensity rises quickly above Jlw ~ 10^^^ erg cm^^ s^-*^ Hz^-*^ sr^^ and the dissociation of H2 in 
the IGM happens abruptly. In this case, the dissociation spheres grow fast enough^ to cover a 
large volume of the IGM before the contribution of distance sources to the dissociating background 
builds up substantially. In contrast, if the dissociating radiation emitted by the sources is small, the 
dissociation spheres grow slowly, while the the additive contribution of distant sources builds up the 
intensity of the dissociating background more quickly. In this case, the dissociation spheres remain 
smaller than the PFRs (and therefore invisible) until the dissociating background has destroyed all 
the H2 in the IGM. 

Figure 7 is analogous to Figure 6, except that wc now show a zoomed region (0.125 x 0.125 
h^^ Mpc) around the most massive object in the 64L05p3 box. This simulation has L^ox = 0.5 

Mpc, = 0.2, (ei7i//47r) = 1.6 x 10~^, {fesc) = Ij and sources with a Population III SED. In 
this sequence of four slices (at z = 17.3, 12.2, 11.3, and 10.2) we recognize the two main processes 



^In Ricotti et al. (2001) wc provide an analytic expression for the cornoving radius, Rd, of the dissociation sphere 
produced by a source that turns on at a = Zi as a function of time: Rd oc {zi + l){SLwty^^- 
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Fig. 5. — Background specific intensity in 10~^^ erg cm~^ s~^ Hz~^ sr~^ units, J2i{i^), as a function 
of photon energy for the 64L05p3 run (Population III) at z = 12. The lower panel shows a zoom of 
the spectrum in the Lyman- Werner bands. In this example, H2 and H Lyman series line opacities 
reduce the intensity of the dissociating background in the Lyman- Werner bands by about one order 
of magnitude (lower line compared to the upper one). 
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Fig. 6. — Slices through the most massive object in the simulation of the 64L05p2 box at z = 19.4 
and 18.5. The box size is L^ox = 0.5 comoving Mpc. Each one of the two panels shows logXHi 
(top-left), logXi^j (top-right), log(l + 6) (bottom-left) and logT (bottom-right). 

that create H2 in the filaments. In the top- left frame at z = 17.3 we can recognize a PFR as a 
shell of H2 surrounding the Hii region that is barely intersected by the slice. In the bottom-left 
frame (z = 11.3) two Hii regions are clearly visible. Inside the Hii regions, the H2 is destroyed. In 
the bottom-right frame {z = 10.2) the Hii regions are recombining (demonstrating that the SF is 
bursting) and new H2 is being reformed inside the relic Hii regions. A more fine inspection*^ of the 
time evolution of this slice shows that at least five H 11 regions form and recombine between z = 20 
and 10 in this small region of the simulation. 

If e*{euv / ^'^){fesc) *C 3.2 X 10~^ the SF mode is less violently bursting, and the radiative 
feedback does not suppress the SFR with respect to the case without radiative transfer. In this 
case, the positive feedback should be dominated by PFRs preceding the ionization fronts rather than 
by PFRs inside relic H 11 regions. From the inspection of some 3D movies showing the formation 
of Hii regions, it appears that galaxy formation could be triggered by the presence of neighboring 
galaxies, in a chain-like process. It is difficult to prove quantitatively this effect, given that galaxies 
tend to be clustered and, in the early universe, galaxies of very different mass can virialize at the 
same redshift. 



^Movies of 2D slices and 3D rendering of the simulations are publicly available on the web at the URL: 
http://casa.colorado.edu/~ricotti/MOVIES.html 
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Fig. 7. — Same as in Figure 6 except for a zoomed region of 0.125^ Mpc^ around the most 
massive object in the 64L05p3 simulation. In this time sequence of images (top: z = 17.3, 12.2 from 
left to right; bottom: z = 11.3, 10.2 from left to right) we recognize the two main processes that 
create H2 in the filaments: PFRs in front of Hll regions and the reformation of H2 inside relic Hll 
regions. The bursting mode of the SF is evident from the continuous formation and recombination 
of the Hii regions in the time sequence of the slices. 
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3.2. Self-Regulated Star Formation History 



In Figure 8 (left) we compare the SFR of three simulations with Lt,ox = 0.5 h ^ Mpc, 
(€[/V'/47r) = 1.6 X 10-^ if esc) = 1, and a Population III SED (64L05p3, 64L05p3b, 64L05p3c), 
when we reduce e* = 0.2 by a factor of 10 and 100. It appears that the SFR is fairly insensitive to 
the value of e^, if we consider a realistic range of values 0.2 < e* < 0.02 . As is reduced from 0.2 to 
0.02, the oscillations of the global SFR with rcdshift become more smooth, but its redshift-averaged 
value gets only a factor of two smaller. When we reduce from 0.02 to 0.002, the global SFR gets 
a factor of 5 smaller. It is important to note that, when the SFR is dominated by "large-halo" 
objects, the global SFR is proportional to e*. This result is based on the comparison of several 
simulations that differ only on the value of e*, some presented here and some presented in Gnedin 



The Hi ionizing background, Jhi, has the same behavior as the SFR (i.e., is insensitive to e*). 
In Figure 8 (right) we show the comoving mean free path^ of Hi ionizing photons Xcom for the 
same simulations shown in Figure 8 (left). It is striking that, even if we change by two orders of 
magnitude, the mean free path of ionizing photons oscillates around a constant value, Xcom^ shown 
in the figure by the solid line given by 



Figure 9 is analogous to Figure 8, but compares simulations varying the value of euv The 
first three runs of the list shown in the figure (thick lines) have the same (fesc) = 1 and e* = 0.2, 
but ejjv is reduced by factors of 10 and 100; the solid, short-dashed and long-dashed lines have 
(et/v'/47r) = 1.6 x 10~^,1.6 x 10^^ and 1.6 x 10~^, respectively. The last three runs of the list 
(thin lines), have e* = 0.02, and varying ejjv] the solid, short-dashed and long-dashcd lines have 
{euv/'^Tr) = 1-6 X 10"^, 1.6 X 10"'^ and 1.6 x 10"^, respectively. Figure 9 (left) shows that the SFR 
is approximatively inversely proportional to euv and insensitive to e*. The inverse proportionality 
relation is evident when the SF is smooth (compare the thick and thin solid lines). When the SF 
is bursting the comparison of different simulations is more difficult but the inverse relation appears 
to hold at least in a limited range of the parameter space. Figure 9 (right) shows that Xcom is 
constrained to not exceed a critical value, A^^^, shown by the solid line of equation (1). Analogous 
to Xcom, Jhi is insensitive to the choice of the free parameters of the simulations. The number of 
ionizing photons that escape in the IGM is propotional to the parameter combination euv {f esc)- 
Therefore changing euv is the same as changing (fesc)- The only difference is that the SED has 
more dissociating photons if we reduce (fesc) insted of euv- In § 3.2.1 we show that, unless (fesc) 



^By definition, Xcom = [c/fcHi(l-|-2)] where km is the mean absorption coefficient weighted by the photoionization 
rate. Neglecting the terms on the left hand side of equation (4) in paper I, we have km = Sm/ Jni- Therefore we can 
derive Xcom from the emissivity Sm (that is proportional to the SFR, {fesc) and euv) and the ionizing background 
intensity Jhi- 



(2000). 




(1) 



DRAFT: February 1, 2008 



19 




Fig. 8. — Global SFR (left) and comoving mean free path of Hi ionizing photons, Xcom (right) as 
a function of redshift for the 64L05p3, 64L05p3b and 64L05p3c runs. The solid, short-dashed, and 
long-dashed lines have = 0.2, 0.02, and 0.002, respectively. The three simulations have L^ox = 0.5 
h-^ Mpc, (euv/^Tr) = 1.6 x W-\ (fesc) = 1, and a Population III SED. The global SFR, Jhi, and 
Xcom, are almost insensitive to for a reasonable range of SFEs (0.2 < e* < 0.02). 

is very small, the dissociating radiation does not affect the SFR. In this regime is the value of the 
parameter combination ejjvifesc) that regulates the SFR. 

If the product e*euv{fesc) is large, Xcom has large oscillations around the critical value, and 
when the product is small, the oscillations are small. Large oscillations of Xcom are associated 
with a strong bursting SF mode. The feedback works in such a way that when Xcom exceeds the 
critical value the SF is suppressed, and consequently Xcom becomes smaller than the critical value. 
It is easy to show that Xcom is related to the the Hii region radii. Therefore, the mechanism that 
self-regulates the SF in "small-halo" objects is related to the size of Hii regions, rather than to the 
intensity of the dissociating background as previously thought. 

Following Gnedin (2000), we can express Xcom as a function of the mean radius, /^hiu of 
the Hii regions well before the overlap phase. By definition, Xcom = [c//c(l + z)] where k is the 
mean absorption coefficient weighted by the photoionization rate F. The mean optical depth for 
Hi ionizing radiation inside an Hii region is (r) = Ruu/ Xcom = (chi^Hi), where uhi is the Hi 
ionization cross section and A^hi is the Hi column density of the halo. The Hi column density of 
a virialized halo is 

Nm = xuirivirRvir ~ 10^^ cm 
Here Uyir is the virial density, Rmr is the virial radius of a M^m Mq halo, and xhi 0.05 
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(derived from the simulations) is the hydrogen neutral fraction. Assuming a frequency-averaged 
value of the Hi photoionization cross section (cthi) ~ 10^^^ cm^, we have: 

Tjcom T>com / of) \ ^ 

= Tf = " ^^'^ (m) ■ 

where -RffJ^ {h~^ kpc ] is the mean free path of ionizing photons. Comparing equation (2) with 
equation (1), we find that the mean radius of the Hii regions produced by "small-halo" objects is 
-^HiT ~ ^ kpc, about the size of the dense filaments and the virial radii of the halos. 

Qualitatively, this result is already evident in Figure 2, which shows that the Hii regions 
remain confined inside the filaments. The PFRs produced ahead of ionization fronts and the relic 
Hii regions continuously reform H2 inside the filaments. The H2 abundance remains high in the 
filaments, even when the dissociating background intensity is sufficiently strong to dissociate all 
the H2 in the lower-density IGM. 
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Fig. 9. — Same as in Figure 8 for the 64L05p2noLW-2, 64L05p2-l, and 64L05p2 runs (thick 
solid, short-dashed, and long-dashed lines) with e* = 0.2, {fesc) = l, and (e[/y/47r) = 1.6 x 
10^'', 1.6 X 10^^, 1.6 X 10^^. The thin solid, short-dashed, and long-dashed lines, show the 64L05p3b- 
1, 64L05p3b, and 64L05p2noLW-hl runs that have e* = 0.02, {fesc) = 1, and (e[/y/47r) = 1.6 x 
10~^, 1.6 X 10~^, 1.6 X 10~^. The three curves in each set of simulations have values of ejjv increased 
by one order of magnitude. It appears that the SFR is about inversely proportional to euv and 
independent of e*. The comoving mGcin frGG pcttti of Hi ionizing photons, Acom? 

and Jhi, instead, 

remain constant. 
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3.2.1. Star Formation History: a Function of {fesc) 

We have seen in the previous section that the SFR is fairly insensitive to over the range 
0.02 < e* < 0.2 and inversely proportional to euvifesc)- The free parameter euv depends mainly on 
the stellar IMF and slightly on the stellar metallicity. Assuming a Salpeter IMF we find {euv / 47r) = 
1.1 X 10"^ for a Population II SED and (e[/y/47r) = 2.5 x 10"^ for a Population III SED. Some 
theoretical arguments suggest that, at high redshift, the IMF could be flatter than a Salpeter IMF. 
If this is the case, we would have (e[/v'/47r) ^ 10~^. The possibility of having a steeper IMF, and 
therefore (e(/v'/4-7r) <C 10^^, is not supported by any theoretical work or observation. 

The value of {fesc) is unknown, and has been an argument of debate form many years. In 
literature, values of (fesc) between 0.5 to zero have been proposed. Theoretical work on (fesc) 
for "small-halo" objects at high-redshift (Ricotti & Shull 2000; Wood & Loeb 2000) finds that 
(fesc) should be very small (10~^ <{fesc) < 10~^) and decreasing with increasing halo masses. 
Observations of (fesc) in nearby starburst galaxies find values of (fesc) ~ 10% (Leitherer et al. 
1995; Hurwitz, Jelinsky, &: Dixon 1997; Heckman et al. 2001; Dcharveng et al. 2001), in agreement 
with theoretical studies (Dove, Shull, Sz Ferrara 2000). Numerical simulations of the reionization 
of the IGM usually adopt a value of (fesc) ~ 1 in order to reionize the universe before redshift 
z ~ 6. A recent study ou Lyman break galaxies (Steidel, Pettini, & Adelberger 2001) agrees with 
numerical simulations in finding {fesc) ~ 0-5 — 1. But there is no observational constraint on the 
value of {fesc) for "small-halo" objects. Reasonably, "small-halo" objects should have smaller {fesc) 
than "large-halo" objects (normal galaxies) since their halos are not collisionally ionized. On the 
other hand, if a substantial fraction of the ISM is photocvaporated or blown away by stellar winds 
and SNe, {fesc) could be larger. We note that {fesc) in this paper is defined as the escape fraction 
of ionizing photons from the resolution element; therefore it is resolution dependent and generally 
larger than {fesc) from the galactic halos. 

In the previous section we showed that, if {euv /^'^){ fesc) ^ 10~^ and e* ^ 0.02, the SFR is 
suppressed by the feedback of ionizing radiation (but not from the dissociating background) . Since 
it is unlikely that (e;7y/47r) < 10~^, the SFR can be increased if {fesc) < 1- As we decrease {fesc), 
the SFR increases almost linearly up to a maximum value determined by the SFR without any 
feedback. This value is proportional to and to the mass resolution of the simulation. In higher 
mass resolution simulations, the number of "small-halo" objects that form stars is larger since we 
resolve many more small-mass objects. We believe that our higher resolution simulation, with 
Lbox = 1 h"^ Mpc, 256^ cells, and mass resolution Mdm = 4.93 x 10^ Mq, is close to fully 
resolving SF for the case without radiative feedback (we need to resolve each halo with about 100 
DM particles). 

Clearly, if {fesc) = 0, there should be no positive feedback, but only the negative feedback of 
the dissociating background, which should determine the SFR. Indeed, if we decrease the value of 
{fesc) below a critical value, the SFR, after reaching the maximum, starts decreasing. This effect is 
shown in Figure 10 (left), where the lines show simulations with = 0.2, Population II SED and 
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ifesc) = 1,0.1,0.01,10-3,10-5 (64L05p2, 64L05p2-lf, 64L05p2-2f, 64L05p2-3f and 64L05p2-5f runs 
respectivelly). At z = 15.5, in the simulation with {fesc) = 0.01, the intensity of the dissociating 
background is about ~ 5 x 10-^^ erg cm-^ Hz-^ sr-^. The value of Jlw is high enough 
to decrease the SFR with respect to the case without radiative transfer by about a factor of two. 
In Figure 10 (right), we show the different importance of the dissociating radiation feedback for a 
Population II and a Population III SED. All the simulations in Figure 10 (right) have {fesc) = 0.01 
(runs 64L05p2-2f, 64L05p3-2f, 64L05p2-2fa and 64L05p3-2fa in Table 1). 
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Fig. 10. — Stellar fraction, /star, as a function of redshift. (Left) The two thin solid lines show jstar 
for a simulation without radiative transfer and without either radiative transfer or H2 cooling. The 
thick lines show simulations with (/esc) = 1, 0.1, 0.01, IQ-^, 10-^ as shown in the label. Positive 
feedback increases the SFR as we decrease (/esc) from 1 to about 10^^. If (/esc) 10-^ the disso- 
ciating background is effective in suppressing the SF. (Right) The thin solid line shows /star for a 
simulation without radiative transfer and e* = 0.2; the thick solid and dotted lines show simulations 
with (/esc) = 0.01, and = 0.2, using a Population II or Population III SED, respectively. The 
thin dashed, thick dashed, and thick dot-dashed lines are analogous to the thin solid, thick solid, 
and thick dotted lines, but for e* = 0.05. 

The effect of decreasing excessively the value of (/esc) on the SFR is twofold: (i) it reduces 
the positive feedback of the EUV radiation; and (ii) it increases the relative importance of the 
dissociating radiation over the ionizing radiation by modifying the SED, g^. The negative feedback 
of the dissociating background starts to affect the SFR, depending on magnitude of the jump, /?, 
at the Lyc frequency of the SED, and the magnitude of the positive feedback. The jump /3 is 
inversely proportional to (/esc) and depends on the SED (intrinsic jump). For the same (/esc), the 
Population II SED produces a jump 10 times larger than the Population III SED. We find that, if 
0.02 ;S ;S 0.2, the critical value for which the negative feedback of the dissociating background 
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starts to affect tfie SFR is {fesc) ~ 10~^ for the Population II SED, and {fesc) ~ 10~^ for the 
Population III SED. 

In Figure 11 we show a case where the positive feedback produces enhanced global SFR with 

respect to the case without radiative feedback. The simulation (thick dashed line) has = 0.02, 
(/g^^) = 1, (et/v'/47r) = 1.6 x 10"^ and a Population III SED (64L05p3b-3 run). The thin solid 
line shows the analogous simulation without radiative transfer. The thick dashed line shows a 
simulation with the same value of the parameter combination euvifesc) as the 64L05p3b-3 run, but 
with Salpeter IMF ((e[/y/47r) = 1.6 x 10"^) and {fesc) = 10~^. The SFR in these two simulations 
should be identical, since it depends on the parameter combination euvifesc)- Instead, the negative 
feedback of the dissociating background reduces the SFR by a factor of three with respect to the 
64L05p3b-3 run. The dissociating background in the simulation with {fesc) = is three orders 
of magnitude higher {Jlw ~ 0.5 x 10~^^ erg cm~^ s~^ Hz~^ sr~^ at z = 11.5) than in the simulation 

with {fesc) = 1- 

The simulations presented in this section have mass resolution Mdm = 3.94 x 10^ h'^ M©. The 
global SFR (without radiative feedback) is underestimated because we do not resolve the lowest 
mass objects that can form stars. Higher resolution simulations could produce an enhanced SFR 
with respect to the case without including any feedback. Indeed, positive feedback can trigger SF in 

halos with a mass smaller than the minimum mass derived in the absence of feedback. Theoretically, 
the lower limit for the DM mass of an object that could form stars is determined by the filtering^ 
mass that is Mp ~ 10^ - 10*^ M© in the redshift range 20 < z < 100. In objects with DM 
masses smaller than Mp, the baryons cannot virialize. 

3.2.2. How does the Self-Regulation Work? 

Summarizing, the feedback prevents the size of H II regions from exceeding R^^ ~ 5 kpc 
(about the size of the dense filaments) . Indeed when the H II regions get bigger than the filaments 
molecular hydrogen is destroyed and the SF is suppressed. On the contrary, when the Hll regions 
are smaller than the filaments or when they recombine after a burst of SF, the dense filaments are 
only partially ionized and the formation rate of molecular hydrogen is maximized. Since galaxy 
formation takes place only in overdense regions, the SF is self-regulated to maximize H2 formation in 
the filaments. Thus, the volume filling factor of the Hll regions remains small. As a result, "small- 
halo" objects cannot reionize the universe. Molecular hydrogen is continuously reformed, in shells 
preceding the Hii regions and in shells inside relic Hll regions, the PFRs found by Ricotti et al. 
(2001). SF is bursting, since it is self-regulated by the two above-mentioned feedback mechanisms. 
The SFR does not depend on e* or on the source SED, but only on {fesc) and the IMF through 



The filtering mass is an averaged Jeans mass that depends on the thermal history of the gas. If the temperature 
remain constant with time, the filtering and Jeans mass are identical. 
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Fig. 11. — Stellar fraction, fstar as a function of redshift for the 64L05noRAD (without radiative 
transfer, thin sohd line), 64L05p3b-3 (thick solid line), and 64L05p3b-3n (thick dashed line) runs. 
The simulations have L^ox = 0.5 Mpc, = 0.02, and a Population III SED. The 64L05p3b-3 
has {fesc) = 10"^ and (e[/y/47r) = 1.6 x 10"^. The 64L05p3b-3n run has (/esc) = 1 and (ei7y/47r) = 
1.6 X 10~^ (flat IMF). The dissociating background is three orders of magnitude higher in the run 

with (fesc) = 10-3. 
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euv If ^uvifesc) is small, the SFR is high, and vice versa. Indeed, if few ionizing photons (per 
each baryon converted into stars) escape from the halo, more stars have to be formed in order to 
produce Hll regions of the size of the filaments. If the product e*€uv{fesc) is large, the SF burst 
is so fast that Hii regions will expand outside the filaments before recombining. This produces a 
temporary halt in the global SF, which appears as a sequence of strong bursts. It is fascinating 
that the SF history of "small-halo" objects depends primarily on a single parameter, €ijv{fesc)- 
This happens because the feedback mechanism acts on a cosmological (rather than galactic) scale. 
We notice, though, that (fesc) should depend slightly on e* (Ricotti & ShuU 2000). 

In Figure 12 we show a sketch of the relative importance of positive and negative feedback from 
EUV (ionizing) and FUV (dissociating) radiation on the SFR. The thick and the thin solid curves 
show the SFR with and without radiative transfer, respectively, for a simulation with = 0.2, 
Population II SED ((€c/v/47r) = 1.1 x 10"^), and (fesc) = 0.01. Without radiative feedback, the 
SFR is proportional to e*. Therefore, for a simulation with e* = 0.02, for instance, the SFR shown 
by the thin solid curve would be a factor of ten smaller. The solid line shows the region where the 
positive feedback from EUV radiation is maximum. The SFR of this line is inversely proportional 
to the parameter combination euv {fesc)- Above the solid line, we have increasing negative feedback 
from EUV radiation (Hll regions become larger than the filaments), and below the line the positive 
feedback from EUV radiation becomes increasingly weak. If the parameter combination e^euvifesc) 
;S 10~^ — 10~^, the SFR is not suppressed by EUV radiative feedback (in this case the solid line lies 
above the thin solid curve). When the flux of ionizing photons in the filaments is too small, either 
because {fesc) is small or because of a low SFR, positive feedback effects are weak. The negative 
feedback of the dissociating background becomes, therefore, important. Below the thin dashed 
line, the negative feedback from FUV radiation starts to dominate the positive feedback from EUV 
radiation. The region between the solid and dashed lines is where the positive feedback dominates 
over negative feedback, and its size is inversely proportional to the jump /3 at the Lyc of the source 
SED. The effect of the FUV negative feedback is usually to delay SF at high-redshift, when the 
SFR, and therefore the ionizing photon flux is lower. At higher redshift, as the number of "small- 
halo" objects increases, the positive feedback dominates and the SFR becomes self-regulated by the 
EUV radiation. In the extreme case of a very small {fesc) {{fesc) ~ 10~^ for a Population II SED, 
and {fesc) ~ 10~^ for a Population III SED), the negative feedback of the dissociating background 
can suppress "small-halo" object formation. 



3.3. Metal Enrichment of the IGM 

Metals are produced by the same massive stars that produce ionizing radiation. Understanding 
the metal enrichment of the low density IGM is a challenging task. Observational constraints from 

the metallicity evolution of the Lya forest allow us to test the models. In Figure 13 we show a 
3D rendering of the IGM metallicity (in solar units) for the 64L05p3 run (L^^^ = 0.5 Mpc, 
e* = 0.2, (€[/i//4vr) = 1.6 x 10~^, {fesc) = 1) and a Population III SED). The opacity and the color 
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Fig. 12. — Sketch showing the relative importance of positive and negative feedback from EUV 
(ionizing) and FUV (dissociating) radiation on the SFR. The thick and thin sohd hnes show the SFR 
with and without radiative feedback, respectively (both for = 0.2 and {fesc) = 0.01). The shaded 
strip is where positive feedback dominates. Above the strip, the SF is suppressed by feedback from 
ionizing (EUV) radiation, and below from the feedback from H2 dissociating radiation. 
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coding are proportional to the logarithm of the metallicity {Z/Zq). The comoving volume filling 
factor of the metal enriched gas increases quickly at high-redshift and slows down as the redshift 
decreases. 




Fig. 13. — 3D rendering of the logarithm of the IGM metallicity (in solar units) for the 64L05p3 
run. The four cubes show the time evolution of the volume filling factor of metal-enriched gas at 
z = 15.7,13.3,12.3, and 10.2. 

In Figure 14 (left) we show the ratio {Z/Zq)v/{Z/Zq)m of the volume- to mass- weighted 
mean metallicities as a function of redshift. This ratio is proportional to the filling factor of 
the metal enriched gas. The solid line shows the 64L05noRAD run {L^ox = 0.5 Mpc, e* = 
0.2, (e;7V'/47r) = 1.6 x 10~^, without radiative transfer). The dotted, dashed, and long-dashed 
lines show simulations with (fesc) = 1 and an increasing value of the parameter combination 
e*(ec/y/4vr)(/e<;c) = 3.2 x 10~^3.2 x lO"'^, and 3.2 x 10"^ We remind the reader that, if the value 
of this parameter combination is high, the global SF, Jhi, and Xcom, are strongly oscillating (bursty 
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SF). The filling factor of metal-enriched gas is larger in the simulations with strongly bursting 
SF. This suggests that photo-evaporation of high-redshift "small-halo" galaxies is an important 
mechanism for transporting metals in the low density IGM (voids). 
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Fig. 14. — (Left) Ratio of the volume- and mass-weighted mean metallicities, {Z / Z(7))v / {Z / Zq) m , 
as a function of redshift for the 64L05noRAD (soUd), 64L05p2-2 (dotted), 64L05p3b (dashed), 
and 64L05p2noLW+l (long-dashed) runs. This ratio is proportional to the volume filling factor of 
the metal enriched gas. The lines, from the bottom to the top, have an increasingly bursting SF 
mode. The SF bursts are stronger when the parameter combination e^euvifesc) is large. (Right) 
Mass-weighted mean metallicity, {Z/Zq)m, as a function of redshift. 

Figure 14 (right) shows the mass-weighted metallicity {Z/Zq)hj for the four simulations in 
Figure 14 (left). It can be easily shown that {Z/Zq)m oc euvfstar, where fstar is the mass fraction in 
stars and euv is the energy in ionizing photons per rest mass energy of H atoms {rriHC^) transformed 
into stars (the number of ionizing photons emitted by a population of stars is proportional to the 
number of heavy elements released in the ISM). The parameter ejjv depends on the IMF and stellar 
metallicity. If the feedback of the EUV radiation regulates the SFR, we have fstar oc (euvifesc))^^ ■ 
In this case, we expect {Z/Zq)m oc{fesc) independent of the stellar metallicity and IMF. We 
conclude that, while the SFR and fstars is inversely proportional to the parameter combination 
euv if esc), the mass of metals produced by "small-halo" depends only on (fesc)- In Figure 14 we 
have used {fesc) = 1- 
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3.4. Realistic Scenarios of Cosmic Evolution during the "Dark Ages" 

In this section we show the evolution of the global properties of the universe in our three most 
realistic simulations: the 128Llp2, 128Llp2-2 and 256Llp3 runs. The mass resolution {Mdm = 
4.93 X 10^ Mo for the 256^ and 3.94 x 10^ h'^ M© for the 128^ runs) and box size, L^ox = 1 

Mpc, of these simulations is sufficiently large to resolve the formation of the first "small-halo" 
and "large-halo" objects. In particular, the 128Llp2-2 and 256Llp3 simulations include the effects 
of secondary electrons and the SED modification caused by assuming realistic values of {fescj- We 
have (/esc) = 0.1, e* = 0.1, and metal-free SED ((e[/y/47r) = 2.5 x 10"^) for the 256Llp3 run. 
For the 128Llp2-2 run we assume {fesc) = 0.01, e* = 0.05 and Z = OMZq Population II SED 
{{euv/^TT) = 1-1 X 10-5). The 128Llp2 run has {fesc) = 1, e* = 0.2, (e[/y/47r) = 1.6 x 10"^ 
and Population II SED. Wc believe that the 256"^ run is very close to the limit of numerical 
convergence. Finally, the formal spatial resolution is 156 pc comoving [B^ = 25 is the parameter 
that regulates the maximum deformation of the Lagrangian mesh: the spatial resolution is ~ 
L/box/i^boxB*)] in the 256Llp3 run, 488 h'^ pc comoving {B^ = 16) in the 128Llp2-2 run, and 
781 pc comoving (B* = 10) in the 128Llp2 run. To give a better idea of the scales resolved by 
the simulations, we remind the reader that the comoving core radius of a just-virialized DM halo 
of mass Mdm is Rc ~ 300 pc^Mdm/W^Mq)^^^ (we have assumed a halo concentration parameter 
c = Ryir/Rc = 10). Wc refer the reader to Paper I for details on the physics included in the code 
and in our convergence studies. 

The thick solid lines in Figure 15 show the comoving SFR and fraction of baryons in stars, 
function of redshift for the 256Llp3 run (left panels), 128Llp2-2 run (center panels) 
and 128Llp2 run (right panels). As a comparison we plot the same runs without radiative transfer 
(dashed lines) and the 64LlnoC run (64^ cells and Ljjox = 1 h"^ Mpc), excluding both radiative 
transfer and H2 cooling (thin solid lines). Therefore, thin solid lines show the contribution to 
the global SFR of "large-halo" objects only, and the dashed lines show the SFR of "small-halo" 
and "large-halo" objects without any feedback effect. The main result shown in this figure is 
that "small-halo" objects are an extremely important (or dominant) fraction of the galaxies at 
least until redshift z ~ 10. Contrary to what is widely believed, their formation is not severely 
suppressed by the dissociating background. We showed in § 3.1 that the dissociating background 
has little influence in determining the SF history. In § 3.2 we demonstrated that the mass fraction 
of "small-halo" objects formed depends, instead, only on the value of (fesc) and the stellar IMF. 

The two panels on the left of Figure 15 show the SFR and fstar for the 256^ cells simulation 
which has {fesc) = 0.1; in this simulation the first stars form at z = 30. At z = 15 the fraction of 
stars in "small-halo" objects is fstar — 2 x 10^'^, about 5 times the mass of stars in "large-halo" 
objects. The radiative feedback has reduced fstar by a factor of 10. Note that fstar for "large-halo" 
objects (thin solid line) and for "small-halo" objects without feedback (thin long-dashed line) scales 
with e*. The fraction of baryons transformed into stars of "small-halo" objects, including radiative 
feedback effects, does not depend on e* as long as it is smaller than the the value without radiative 
transfer. Since, for "large-halo" objects, fstar e*, the relative importance of "small-halo" to 
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Fig. 15. — The thick solid hnes show the SFR (top panels) and fraction of baryons in stars, fstar 
(bottom panels), as a function of redshift for the 256Llp3 (left panels), 128Llp2-2 (center panels) 
and 128Llp2 (right panels) simulations. The box size is 1 h^^ Mpc comoving. The 256Llp3 run 
has a Population III SED, (fesc) = 0.1 and = 0.1; the 128Llp2-2 run has Population II SED, 
ifesc) = 0.01 and = 0.05, and the 128Llp2 run has Population II SED, (fesc) = 1 and e,, = 0.2. 
The thin dashed lines show the same simulations, without radiative transfer, and the thin solid lines 
excluding both radiative transfer and H2 cooling (i.e., the contribution of "large-halo" objects). 
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"large-halo" objects is inversely proportional to (which here is 0.1). For instance, if we had 
chosen e* < 0.1 for this simulation, fstar of "small-halo" objects would have been as large as 50 
times fstar of "large-halo" objects at z = 15, and the radiative feedback would have suppressed the 
SF by only a small factor. 

We show an example of such a case in the two central panels of Figure 15. In this 128^ 
simulation we use {fesc) = 0.01 and e* = 0.05. Here, the first stars form at z = 26 because of the 
lower mass resolution. Radiative feedback suppresses the SFR only by a factor of 2. At 2; = 15, 
fstar — 3 X 10~^, and the contribution of "small-halo" objects to star formation is about 10 times 
the contribution of "large-halo" objects. At z = 10, "small-halo" objects still dominate fstar by a 
factor of 3. Note that in this simulation the dissociating background intensity is very high and, as 
shown by the short-dashed lines in Figure 16, the H~ and H2 abundances in the IGM are extremely 
small (~ 10^-*^^ relative to H). 

The two panels on the right show the SFR and fstar for the 128^ simulation, which has (fesc) 
= 1 and e* = 0.2. At z = 15, the fraction of stars in "small-halo" objects is fstar — 2 x 10~^, about 
1/5 the mass of stars in "large-halo" objects. The radiative feedback has reduced fstar by a factor 
of 100. 

In Figure 16 we show the mass- and volume- weighted mean molecular abundances (xj/j); 
{Xfj+ ) and {xff- ) as a function of redshift for the three simulations with radiative transfer shown 
in Figure 15. The volume- weighted abundances are about two orders of magnitude smaller than 
the mass-weighted abundance. This shows that H2 is much more abundant in dense regions. The 
dissociating background destroys the H2 in the low-density IGM in the redshift range 20 < 2; < 
25, depending on the choice of the free parameters in the simulation. In the dense regions, the 
production/destruction of H2 sets its abundance to a quasi-constant value that depends on the 
SED of the sources. At redshift z ~ 15, the abundance of H~, the main catalyst for H2 production, 
reaches its maximum value and starts to decrease. Consequently, the H2 abundance also decreases 
after z ~ 15. 

Figure 17 is analogous to Figure 16, but here we show the mass- and volume- weighted ion- 
ized H and He mean abundances {xnn), {xucn), (a^Hem) (top), and the metallicity (Z/Zq) and 
temperature (T) (bottom). If reionization takes place at redshift z ~ 6, as recent high-redshift 
quasar observations seem to suggest (Becker et al. 2001; Djorgovski et al. 2001), the run shown 
by the long-dashed line reionizes the universe too early, while the run shown by the short-dashed 
line reionizes the universe too late. This observation does not constrain the value of {fesc) or 
for "small-halo" objects, since we expect that {fesc), at least, will vary as a function of redshift. In 
this paper, we assume for simplicity that {fesc), e*, euv, and are constants. 

The mass-weighted metallicity, {Z/Zq)m, is proportional to the total mass of metals produced 

by the stars, and therefore is proportional to the SFR. Since we have shown that "small-halo" 
objects dominate the SFR if {fesc) < Ij at least before redshift z ~ 10, their metal production 
dominates by the same amount. Moreover, we have shown in § 3.3 that, if SF is bursting, photo- 
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Fig. 16. — Mass- and volume- weighted molecular abundances (x/fj); {^h^) ^'^d {x^-) for the three 
runs from Figure 15. The solid hues show the 256Llp3 run {{fesc) = 0.1 and = 0.1), the short- 
dashed hues show the 128Llp2-2 run {{fesc) = 0.01 and e* = 0.05), and the long-dashed lines show 
the 128Llp2 run {{fesc) = 1 and = 0.2). The drops in H2 abundance, while and H~ are stiU 
high, are due to the rising flux of the dissociating radiation. 
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Fig. 17. — (Top) Mass- and volume-weighted ionized H and He abundances: (xhh), (a;Hcn)i (a^Hcm)- 
(Bottom) Mass- and volume- weighted metallicity {Z/Zq) and temperature (T). The solid lines show 
the 256Llp3 run {{fesc) = 0.1 and e* = 0.1), the short-dashed lines the 128Llp2-2 run {{fesc) = 0.01 
and e* = 0.05), and the long-dashed lines show the 128Llp2 run {{fesc) = 1 and e* = 0.2). 
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evaporation of small galaxies is the main process that pollutes the low-density IGM (voids) with 
metals. Finally, SNe explosions, not yet included in the simulations, can transport metals into the 
voids efficiently, since "small-halo" objects are numerous and the voids are small at high redshift. 

In Figure 18 (left) we show the evolution of the Hi ionizing background Jhi (thick lines) and 
the He ii ionizing background Jhcii (thin lines) for the same three simulations shown in the previous 
figure. In Figure 18 (right) we show the comoving mean free path of the Hi ionizing photons Ahi 
(thick lines) and of the Heii ionizing photons Aneii (thin lines). Reionization, defined as the overlap 
of Hii regions, occurs when Ahi ~ Ds/2 0.1 Mpc, where Dg is the mean comoving distance 
between the ionizing sources (Gnedin 2000). 

a a 
0.04 0.06 0.08 0.1 0,04 0.06 0.08 0.1 
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Fig. 18. — (Left) Evolution of the Hi ionizing background Jhi (thick lines) and the Heii ionizing 
background Jhcii (thin lines). (Right) Comoving mean free path of the Hi ionizing photons Am 
(thick lines) and of the Heii ionizing photons Aneii (thin lines). The solid lines show the 256Llp3 
run {{fesc) = 0.1 and e* = 0.1), the short-dashed lines the 128Llp2-2 run {{fesc) = 0.01 and 
= 0.05), and the long-dashed lines show the 128Llp2 run {{fesc) = 1 and e* = 0.2). 



4. Discussion and Summary 

In this paper we studied the the formation and evolution of the first galaxies in a AC DM 
cosmology. Our results are based on 3D cosmological simulations that include, for the first time, a 
self-consistent treatment of radiative transfer. The simulations include continuum radiative transfer 
using the "Optically Thin Variable Eddington Tensor" (OTVET) approximation and line-radiative 
transfer in the H2 Lyman- Werner bands of the background radiation. Chemical and thermal pro- 
cesses are treated in detail, particularly the ones relevant for H2 formation and destruction. The 
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details about the numerical methods and the physics included in the simulations are treated in 

a companion paper (Paper I). In Paper I we have also performed careful convergence analysis. 
It appears that we are very close to the convergence limit when we use a 256^ box. In smaller 
simulations the SF is underestimated. 

The main result is that positive feedback processes dominate negative feedback of the disso- 
ciating background, and therefore star formation in "small-halo" objects is not suppressed. The 
main parameter that determine the importance of "small-halo" objects is {fesc)- If ifesc) ~ 0.01, 
"small-halo" objects dominate the galaxy mass function at least until redshift z ~ 9, and feedback 
produces a bursting SF in these objects. Because the ionization fronts are confined to the dense 
filaments, reionization of the IGM cannot be produced by "small-halo" objects. However, we find 
that "small-halo" objects are important in enriching the low-density IGM with heavy elements. 

The main processes that are not included in our simulations are H2 self-shielding and mechan- 
ical and thermal energy injection caused by SN explosions. H2 self-shielding, depending on the 
choice of the free parameters, could be relevant and should reduce the effects of the dissociating 
radiation. Since we have found that the dissociating radiation has a negligible role in regulating 
SF, we expect that the results of our simulations should not be modified by H2 self-shielding. The 
effects of dissociating radiation become important when {fesc) ~ 10~^. H2 self-shielding effects 
probably reduce the value of (fesc) for which dissociating radiation becomes important to a smaller 
value. Although we have not included them in our simulations, SN explosions can be important; 
their feedback could be responsible for a self-regulating global SF and contribute to spread metals 
into the low-density IGM. Unfortunately, the dynamical and thermal effects of SN explosions on the 
ISM of galaxies and IGM are not fully understood. The investigation of the effects of mechanical 
feedback from SNe explosions will be the subject of our future work. 

We cannot rule out alternative cosmologies, in which "small-halo" objects do not form. For 
example, in warm dark matter (WDM) cosmologies, the free-streaming of the DM particles is large 
enough to suppress the formation of dwarf galaxies at high redshift. It is also possible that the 
initial power spectrum does not have enough power on small scales to form "small-halo" objects. 
Theoretical arguments based on observations of the Lya forest, quasars at high redshift, and stellar 
populations in dwarf galaxies pose some constraints on the mass of WDM particle candidates 
{^wdm ~ 1 keV). Current observations do not rule out these alternative cosmological models; but 
they may not solve the problems faced by CDM. 

In the following list we summarize the results of this work in more detail: 

(1) SFR in "small- halo" objects is self regulated and depends only on {fescjt for fixed ejjv 
(known IMF). The SFR is almost independent of the star formation efficiency, e*, the dissociating 
background, and the SED (metal-free or metal-rich objects) of the sources. It depends only on 

{fesc) and the IMF (number of ionizing photons per baryon converted into stars). If {fesc) is small, 
the SFR is high, while if {fesci 0.01 the SF is only slightly suppressed by radiative feedback. In 
this case, the maximum SFR is proportional to e*, and "small-halo" objects dominate the galaxy 
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mass function at least until redshift z ~ 10. 

(2) "Small-halo" star formation is intrinsically "bursting". Star formation is regulated by 
competing negative and positive feedback from EUV (ionizing) radiation. The H II regions produced 
by "small-halo" objects are confined to the dense IGM filaments. Reionization of the IGM cannot 
happen until massive galaxies are formed. In contrast to massive objects, which reionize voids first, 
"small-halo" objects partially ionize only the dense filaments while leaving the voids neutral. 

(3) Galaxy formation is triggered by the presence of neighboring galaxies. The SFR does 
not depend strongly on the dissociating background intensity or on the H2 abundance in the IGM. 
The self-regulation of star formation relies on H2 being continuously reformed in positive feedback 
regions. The H2 formation happens only in the filaments, both inside relic H II regions that never 
expand into the low density IGM and in the H2 shells just in front of the H II regions (Ricotti et al. 
2001). 

(4) The dissociating (FUV) radiation reduces the SFR only if (fesc) is very small. If 

ifesc) ~ 10"-'^ for a Population II SED or {fesc) ~ 10"^ for a Population III SED, the negative 
feedback of the dissociating background suppresses "small-halo" object formation, more efficiently 
at high-redshift. 

(5) "Small-halo" objects dominate the metal pollution of the low-density IGM. The 

transport of metals from the galaxies to the low-density IGM happens because of the continuous 
formation and photo-evaporation of small mass ("small-halo") halos. The metal production, if the 
SF is self-regulated by EUV radiation, is independent of the SEE, SED, and IMF, but depends 
only on {fesc}- 

In conclusion, we have shown that, if SN explosions do not suppress the formation of "small- 
halo" objects and CDM cosmogonies prove to be correct, "small-halo" objects should have profound 
effects on cosmic evolution. Observations of dwarf spheroidal galaxies, metallicity of the Lja forest, 
and stellar populations in the halo of the Milky Way could verify this model. Computational 
limitations prevent us from evolving a representative sample of the universe to redshifts z < 9. 
At lower redshifts, the bulk of "small-halo" objects merge, forming larger mass galaxies, but some 
of them might survive almost unaffected by the environment. Reionization is probably affecting 
the ISM of these small galaxies quite substantially, photoevaporating the remaining unshielded 
gas. These "fossil" "small-halo" objects could be identified with at least some dwarf spheroidal 
galaxies in the Local Group. In a paper currently in preparation, we study the properties of the 
simulated "small-halo" objects in order to understand whether this link is real. The prediction of a 
large population of "small-halo" objects offers a challenging test to verify CDM cosmogonies. We 
speculate that, in the near future, given the rapid progress of computational power and the rapid 
growth of observational data of cosmological interest, detailed cosmological simulations will allow 
us to constrain the evolution of our free-parameters: {fesc), ^UVi and e*, or possibly the nature of 
the dark matter. 
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